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^ ■ 

■<sj- ■ Abstract 

o : 

^ ■ We describe the observing simulation software FISVI (FIS Virtual Instrument), 

^ ' which was developed for the Far-Infrared Surveyor (FIS) that will be on the Japanese 

Qh| infrared astronomy mission ASTRO-F. The FISVI has two purposes: one is to check 

O I the specifications and performances of the ASTRO-F/FIS as a whole; the other is to 

^ ', prepare input data sets for the data analysis softwares prior to launch. In the FISVI, 

, special care was taken by introducing the "Compiled PSF (Point Spread Function)" to 

•'-^ ■ optimise inevitable, but time-consuming, convolution processes. With the Compiled 

^ ■ PSF, we reduce the computation time by an order of magnitude. The photon and 

readout noises are included in the simulations. We estimate the detection limits for 
point sources from the simulation of virtual patches of the sky mostly consisting of 
distant galaxies. We studied the importance of source confusion for simple power-law 
models for N{> S), the number of sources brighter than 5*. We found that source 
confusion plays a dominant role in the detection limits only for models with rapid 
luminosity evolution for the galaxy counts, the evolution of which is suggested by 
recent observations. 

Key words: methods: data analysis — techniques: image processing — galaxies: 
photometry — infrared: galaxies 
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Table 1. Specifications of the FIS. 





Wavelength range 


Array size 


Pixel size 


Pitch size 


Sampling rate 


Band 


(fim) 


(pixel) 


(arcsec) 


(arcsec) 


(Hz) 


WIDE-L 


110 - 200 


15 X 3 


44.2 


49.1 


15.2 


N170 


150 - 200 


15 X 2 


44.2 


49.1 


15.2 


WIDE-S 


50 - 110 


20 X 3 


26.8 


29.5 


22.8 


N60 


50 - 75 


20 X 2 


26.8 


29.5 


22.8 



1. Introduction 

The FIS (Far-Infrared Surveyor) is one of the focal plane instruments of the ASTRO- 
F mission (previously known as IRIS) (Murakami 1998; Shibai 2000; Nakagawa 2001). The 
ASTRO-F satellite will be launched into a sun- synchronous orbit at an altitude of 750 km, 
which corresponds to an orbital period of 100 min. The telescope, which is cooled down to 
5.1-5.8 K, has a 67 cm primary mirror. The major task of this mission is to carry out an 
all-sky survey across the 50-200 /im range. The basic parameters of the ASTRO-F/FIS are 
summarized in table 1 (see also Kawada 2000). 

ASTRO-F/FIS will bring data with much higher sensitivity and angular resolution than 
those of IRAS (see Kawada 2000 for detailed comparison). Such data sets will be of great value 
for many areas of astrophysics, including cosmology, galaxy evolution, interstellar medium, and 
asteroids. 

Generally speaking, the hardware characteristics of each component in a space mission 
can be measured in the laboratory. However, it is very difficult to make end-to-end tests of 
a mission in the laboratory. Hence, based on data measured for each component, numerical 
simulations are frequently used to understand the instrument performances as a whole (e.g., 
Garcia et al. 1998; Boggs, Jean 2001). Moreover, the complicated interplay between the celestial 
sources and hardware specifications can be studied only by a simulation prior to the launch. 

We have constructed a software simulator called the FISVI representing Virtual 
Instrument of the FIS, that can simulate the data stream of ASTRO-F/FIS (Jeong et al. 2000). 
This work is an extension of initial work by Matsuura et al. (2001). The purposes of the FISVI 
are : (1) to confirm the performance of the hardware as a whole and (2) to generate simulated 
FIS survey data sets as inputs for data-reduction software prior to launch. 

One of the key questions regarding the performance of ASTRO-F is the effective detection 
limit for faint sources. Depending on the size of the sources compared to the beam size of 
ASTRO-F, the source can be either extended or point-like, and the detection limits depend on 
the nature of the sources. In the present work, we only consider point sources. 

There are several factors contributing to the detection limits. The sensitivity of the 
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detectors and the entire telescope system allows only sources brighter than a certain threshold 
to be reliably measured. Since the photons follow Poisson statistics, the background photons 
due to the sky brightness as well as the telescope emission should fluctuate, and a meaningful 
detection of a source can be made only if the signal from the source exceeds the level of the 
fluctuations. The sky confusion noise by the cirrus emission causes an uncertainty in the 
determination of the source flux, due to the variation of the sky brightness (Herbstmeier et 
al. 1998; Kiss et al. 2001). The readout process also adds more fluctuations. Moreover, the 
measurement of the brightness of a source can be further influenced by neighboring sources 
if more than one source lies within a single beam of the telescope. The final detection limit 
should thus depend on the performance of the entire system, the brightness of sky and telescope 
emission, readout process, and the distribution of sources as a function of the flux. 

There have been a number of estimates of detection limits based on the available lab- 
oratory data (e.g., Kawada 1998, 2000) using simple calculations. Clearly, a more realistic 
estimation can be made by using numerical simulations. In the present work, we carried out 
simulations of the ASTRO-F/FIS observations under several different circumstances in order 
to obtain still more reliable detection limits which can be used to design scientific projects. 

The present paper is organised as follows. In section 2, we briefly describe the design 
and the structure of the FISVI. In section 3, we explain how we obtain the observed images 
based on the simulated data set. In section 4, we make estimate on the detection limits of the 
ASTRO-F/FIS under various circumstances. First, we estimate the detection limits of a single 
isolated point source while considering only photon and readout noises. Also, we estimate the 
confusion noise (Condon 1974; Franceschini et al. 1989) for distributed sources using a simple 
formula. By carrying out aperture photometry to the simulated images, we finally obtain 
combined detection limits that include photon, readout, and confusion noises. The final section 
summarises our conclusions. 

2. Structure of FISVI Software 

The algorithm of the FISVI software is shown in figure 1. The input data file provides 
the coordinates and fluxes of the sources in the sky. Although the sources would appear either 
point-like or extended, we concentrate on point sources in this paper. The software first makes 
images on the focal plane by convolving the point sources and the Point Spread Function (PSF) 
of the telescope and the instrument. The software generates time-series data for each pixel by 
simulating the scanning procedure of the ASTRO-F/FIS survey mode observations. 

Since the PSF, the filter transmission, and the detector response depend on the wave- 
length of incoming photons, we need to do repeated calculations (procedures boxed in the left 
panel of figure 1) for different wavelengths within the individual FIS bands, as shown in the 
left panel of figure 1. To elude this and speed up the procedure, we introduce the Compiled 
PSF in this work, with which we can perform this scanning procedure at once, as shown in the 
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FigureOl .ps 



Fig. 1. Flow charts of FISVI. The left chart shows a straightforward procedure based on the realistic 
photon path, where repeated calculations would be necessary over the wavelength grids ("A Convolution"). 
The right chart, on the other hand, shows the accelerated algorithm using the Compiled PSF for the FISVI. 

simplified flow chart in the right panel of flgure 1. A more detailed discussion on the gains in 
the computational time and possible errors due to the use of the Compiled PSF are presented 
in appendix 1. 

The readout values for each pixel are represented by a series of integrated charges taken 
over the area covered by the pixel, sampled at regular time intervals. The integrated charges are 
set to zero at every reset interval. The time series of the integrated charges are differentiated 
to obtain the charges accumulated during the sampling interval (see appendices 2 and 3 for 
detailed process). We also generate the photon and readout noise and include them to be part 
of readout values. A more detailed discussion on the implementation of noise is presented in 
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subsection 4.1. The time-series data are converted into the brightness distribution on the sky, 
and are used to reconstruct the images, as described in section 3. 

3. Image Reconstruction 

The FISVI generates time-series data for each pixel. In figure 2, we show a series of 
readout values of a pixel that scans across a point source. No reset was applied during the 
readout sequence shown in this figure because the reset time interval is usually much longer 
than the passage of a Compiled PSF over a point source. The differentiation (subtraction of 
adjacent sampling points) of this curve gives the signal obtained during a sampling interval by 
one detector pixel, which is shown in the lower panel of figure 2. 

The pixel readouts can be used to reconstruct the images. In the current implementation 
of the FISVI, following method was used to generate the image. In order to reconstruct the 
image, we assume that pixel value represents the uniform intensity over the pixel surface. This 
means that a particular point can be covered by more than one readout. We always take the 
average values of multiple readouts in order to construct images (see figure 3). Due to the 
convolution of the image with the pixel size, the output image will be blurred slightly. 

4. Estimations of FIS Performance 

An estimation of the detection limits for the planned mission is very important. For 
ASTRO-F, the detection limits were estimated by using analytic methods (Kawada 1998, 2000). 
In the present work, we made a numerical estimate for a single point source using the latest 
information for the detectors and filters, and compared them with the photometric results on 
the FISVI generated images that contain a large number of point sources. 

4.I. Detection Limits for a Single Point Source 

The detection limits for a single point source depend on the level of noise. There are 
several sources of noise: photon noise due to the sky background and thermal emission from the 
telescope, and readout noise. The sky background varies significantly from place to place in the 
sky. On average, the infrared sky becomes brighter in the Galactic plane, and diminishes toward 
the Galactic poles. Within the Galactic plane, the emission from the Galactic center direction 
appears to be brighter than towards the anti-center direction. Because of thermal emission by 
interplanetary dust particles, the ecliptic plane is also brighter than the ecliptic pole region. 
In figure 4, we show the assumed surface brightness distribution of background emissions from 
the interstellar dust, the interplanetary dust and the telescope, for the purpose of generating 
photon noises. These background emissions from the sky are assumed to correspond to the 
dark part of the sky and the sky confusion noise due to the structure of the cirrus emission is 
not considered. The telescope temperature is assumed to be 6 K, as a conservative number. 
In figure 4, we also plotted the thermal emission from the 6.5 K telescope as a comparison. 
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Figure02a. ps 



Figure02b . ps 



Fig. 2. Example of a scries of readout values, which corresponds to the integrated charges since the 
last reset [see equation (A5)], of a WIDE-L pixel that scans through a point source (left panel). The 
differentiation of the integrated charges as shown in the right panel corresponds to the signal obtained 
during a sampling interval by one detector pixel passing the image of a point source. The sampling interval 
was 14."2, corresponding to the 15.2 Hz readout (see table 1). 
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FigureOS .ps 



Fig. 3. Schematic figure for image reconstruction by pixel averages. At any given point, we take the 
average of the pixel readouts that were covered by those pixels. In the figure, the darker area means the 
area that was covered more. 

Evidently, the contribution from the telescope is smaller than that from the interplanetary 
or interstellar dust as long as the telescope temperature is lower than 6.5 K for the entire 
FIS bands. The sky brightness throughout the spectral region of the FIS varies from 5 to 7 
MJysr"^. Obviously, we would need to apply a position-dependent background brightness for 
more realistic sky simulations, which affects the photon noise. The incoming photon stream on 
pixels due to background emission is assumed to follow Poisson statistics. 

The readout circuit also generates uncertainties of the output values, called readout 
noise. This type of noise is independent of the sampling rate and the integration time, and we 
assumed the total noise in the effective bandwidth at the first stage of the field effect transistor 
(FET) gate to be 3 /iV. In the simulation, we assumed that the readout noise follows Gaussian 
statistics. 
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Figure04.ps 



Fig. 4. Assumed background emissions. We consider three components for the background emission, i.e., 
interstellar dust (dotted), interplanetary dust (dashed) and telescope emission assuming 6 K (dot-dashed 
line) or 6.5 K black body (long dot-dashed line). In our simulations, the telescope temperature is always 
assumed to be 6 K. 
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Table 2. Simple estimates of 5cr detection of single pixels and the ratios of photon-to-readout noises. 





5cr Dot. Limit (mJy)* 




WIDE-L 


39 


1.3 


N170 


76 


1.8 


WIDE-S 


20 


1.6 


N60 


52 


2.5 



* Average flux density in the bandwidth. 
^ Readout-to-photon noise ratio. 

4.1.1. Simple estimation 

The sky brightness throughout the spectral region of the FIS varies from 5 to 7 MJysr"^. 
The integrated photons fluctuate following the Poisson statistics while the readout process adds 
readout noise, which is assumed to follow Gaussian statistics. The r.m.s. fluctuation of voltage 
across the integrating charge due to readout noise can be converted to the fluctuation in the 
number of charges by 

V 

-"-^rms ^ 1 V y 

where C is the capacitance of the charge integrators [7 pF for SW (short wavelength) and 10 
pF for LW (long wavelength) bands, respectively], and e is the elementary charge. The total 
noise is a combination of photon and readout noise. 

If we assume that a single pixel detector receives the entire photon flux of the point 
source, we can obtain the accumulated charge during 'the effective integration time' that elapses 
until the detector pixel passes through one point. For a photoconductor, the noise by this 
photon flux arises from the sequence of generations and recombinations of photoelectrons. We 
calculated this generation-recombination noise (G-R noise), /g-r (Rieke 1994) using 

= Ae'^vG'df, (2) 

where (f is the photon flux, t] is the quantum efficiency, G is the photoconductive gain, and df 
is the effective bandwidth. We assumed that the source has the SED of a 40 K blackbody. The 
5a detection limits computed in this way for all FIS bands are shown in table 2. Also shown in 
this table is the relative importance of the photon and readout noise. In all cases, the readout 
noise is more important than the photon noise, with narrow bands (N170 and N60) being more 
dominated by readout noise. 

4.1.2. Estimation using scanning simulations of a single pixel 

We also estimate the detection limits from the detector scanning routines in the FISVI 
for a single pixel. The behavior of the readout values as a function of the sampling sequence 
is shown in figure 2. The contribution due to background can be obtained by subtracting the 
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Table 3. 5cr detection limits of FIS bands from scanning simulations with a single pixel. 



Detection Limits (mJy) 



WIDE-L 


40 


N170 


80 


WIDE-S 


20 


N60 


47 



contribution from the source alone. The expected amount of the fluctuation is proportional to 
Gy/ipr] for a given span of the scanning period of ti to ^2 [see equation (2)]. The total amount 
of fluctuation of the readout value due to noises during the same scanning span, atot, is 

f^tot = V^f^ph + ^r. (3) 

where Uph and are the fluctuation due to the photon and readout noise, respectively. Here, we 
assume that the readout noise is always a constant while the amount of charge fluctuation due to 
the photon noise increases as Gy/(fr], as dictated by the Poisson nature. For a given brightness 
of a source, we can obtain the S/N ratio if we specify ti and ^2- Since the signal (photocurrent) 
and the photon noise are proportional to G, S/N depends on y/fj on the condition that the 
photon noise is the dominant case. From equation (A7) and the assumption G = 0.9, we can 
obtain the quantum efficiency, 77, as 0.17 for SW and 0.27 for LW detectors, respectively. The 
determination of ^2 and ti was done to maximize the S/N. We find that this can be done when 
we start the scanning at a distance of 2VFh and continue until the same distance in the opposite 
side, where Wn is the full width at half maximum of the beam patterns (see subsubsection 
4.2.2 for details). The 5a detection limits determined in this way for all FIS bands are listed in 
table 3. These estimates also assume a blackbody source with a temperature of 40 K. We find 
that the estimates using the simple method described in subsubsection 4.1.1 and here agree 
very well each other. The largest discrepancy occurs for the N60 band, where the estimated 
detection limit using scanning simulation is lower by around 10%. The instrumental noise 
in ISO observation is estimated to be 15-45 mJy (Herbstmeier et al. 1998; Dole et al. 2001). 
Assuming our background brightness of ~ 5 MJy sr~^, this noise level is similar with our 
estimation in the wide bands. We analyse the photometric accuracy of point sources in more 
realistic simulations with distributed sources below. 

4-2. Simulations with Distributed Point Sources and Realistic Detector Configurations 

The FISVI takes into account the full configuration of FIS detector arrays. We now 
discuss the simulations over a finite patch of the sky with randomly distributed sources. By 
carrying out the photometry of simulated images, we should be able to determine more realistic 
detection limits. 
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Most faint sources to be observed by the ASTRO-F/FIS are expected to be distant 
galaxies. Since the size of the PSFs at far-infrared wavelengths is relatively large, we expect 
that the number of sources overlapped within a given PSF will be larger. In such a situation, 
the source confusion would be important for faint sources. In this section, we consider how the 
source confusion would affect the observations by the ASTRO-F/FIS. 

4-2.1. Source distribution 

The effect of confusion depends on the distribution of sources in the sky and the PSF. 
We assume that A^(> S), the number of sources whose flux is greater than flux S, as a power-law 
on S, 



for S^in < S < Smax, where A^o and 5*0 are normalisation constants. For uniformly distributed 
sources in Euclidean space, 7 is 1.5. If the galaxies experience strong luminosity evolution from 
active to less active star formation with time, 7 will become greater than 1.5. The curved space 
could also give 7 different from 1.5. The analysis of IR galaxy counts by ISO and SCUBA 
suggests that 7 would be greater than 1.5 but lower than 2.5 at around ~ 150 mJy (Puget et 
al. 1999; Franceschini et al. 2001; Pearson 2001; Dole et al. 2001). Matsuhara et al. (2000) 
suggested that 7 could be steeper than 2.5 based on the fluctuation analysis due to the strong 
evolution. In this paper, we examine three cases: 7 = 1.5, 2.5, and 3.0. We fixed Smin = 10 mJy 
throughout the paper. Since there is no divergence due to Smax, "we do not fix this number. 

We need to specify the normalisation constants, A^q, at a given flux Sq, which is set to 
be 100 mJy. These constants are determined from IR galaxy counts normalised to Euclidean 
law [A^(> S) oc S"^'^] at 90 fim based on the IRAS survey and the European Large Area ISO 
Survey (Efstathiou et al. 2000; Franceschini et al. 2001). In the following cases, though the 
source count results are different for different bands and galaxy evolution, we assumed that there 
are 10 sources brighter than 100 mJy per square degree, i.e., A^o(> 100 mJy) = 10, in every 
observational band and the SED of all sources are fiat within a given FIS band. The number 
density of sources was estimated to be 316 per square degrees corresponding to 0.2 within a circle 
of radius of Wn in LW bands for 7 = 1.5 with the above normalisation. The density becomes 
10-times larger for the case of 7 = 2.5 and the case of 7 = 1.5 and Nq = 100, and 19 times larger 
for the case of 7 = 3.0 and A^o = 60. We expect that source confusion becomes important for 
these distribution. The distribution of sources in the sky is assumed to be uniform Poisson. 
In this work, we want to check the pure confusion effect for the same distributed galaxies by 
excluding other factors, e.g., various types of SED, the redshift distribution, the luminosity 
function and the galaxy evolution. For a comparison, we also check other cases: the Euclidean 
space with a large normalisation constant (Aq = 100) and an extreme case ( 7 = 3.0, Nq = 60) 
(Matsuhara et al. 2000). 




(4) 
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4-2.2. Simple estimate of the confusion noise 

Although the clustering of sources could also affect the confusion noise, we ignore such 
a possibility for simplicity. Following Condon (1974) and Franceschini (1989), we obtain the 
noise due to confusion as 



where x [= Sh{6,(j))] is the intensity, Xc is a cutoff value, and R{x) is the mean number of 
sources within the normalised beam pattern, h{6,(j)): 



where ^(5*) is a differential number count. 

In this calculation, we use the beam pattern (see figure 5), which is obtained from a 
simulated image of an isolated point source using the FISVI without noises. The beam pattern 
obtained in this way is somewhat wider than the Compiled PSF due to pixel convolution. We 
also use the differential number count obtained from the same source distribution assumed in 
subsubsection 4.2.1. These considerations are for the purpose of comparing with the results 
from the photometry in subsection 4.3. We list the 5a confusion noise in table 4, obtained by 
using equation (5) for 7 = 1.5, 7 = 2.5, and 7 = 3.0. We also estimated the crowded fields for 
7 = 1.5 by simply increasing Nq by a large factor, i.e., A'^o(> 100 mJy) = 100. The 5a confusion 
noise is the same for the wide and narrow bands, because the beam patterns are similar for two 
bands. Because of differences in the size of beam profiles between long and short wavelengths, 
the detection limits for LW are higher than those of SW bands. The detection limit by confusion 
is approximately proportional to A'^^^'*^. 

The confusion noise in FIRBACK survey by ISO is estimated to be around ac — 45 
mJy (Dole et al. 2001). In our case, we used the slope of the source distribution as 7 = 1.5 
or 7 = 2.5 and set the normalisation constant as A^o(> 100 mJy) = 10 by using the 90 fim 
source count result (Efstathiou et al. 2000). Though the slope of the source count by Dole 
et al. (2001) is similar to the Euclidean space (7 = 1.5), the normalisation constant should 
be different because the source density and the galaxy evolution is different in other bands. 
Therefore, these discrepancies result from the different normalisation and the cutoff flux (S'^in 
= 10 mJy). 

4.3. Realistic Simulations 

The assumed source distribution of equation (4) can be used to simulate the observed 
sky by the ASTRO-F/FIS. By analysing the simulated images, we can address the effects of 
the various sources of noises to the observation in a more realistic way. 




(5) 




(6) 
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Figure05.ps 



Fig. 5. Comparison between the Compiled PSF (dotted) and the beam pattern (sohd hne) used in calcu- 
lating the theoretical confusion (WIDE-L). Because we assumed a flat SED for all sources in this simulation, 
we used one Compiled PSF in the PSF- convolution. 
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Table 4. Scr detection limits due to confusion noise based on theoretical estimates. 



Band 


7= 1.5 

iVo* = 10 
(mJy) 


7 = 1.5 

iVo* = 100 
(mJy) 


7 = 2.5 
NS = 10 
(mJy) 


7 = 3.0 
= 60 
(mJy) 


WIDE-L 


23 


108 


50 


196 


N170 


24 


115 


52 


204 


WIDE-S 


12 


54 


35 


123 


N60 


11 


52 


34 


121 



* No(> 100 mJy). Number per square degree. 

4-3.1. Realistic simulations 

Using the FISVI, we generated two-dimensional images in the FIS bands for two different 
virtual sky data with different 7. We made two different sets of simulations. One was with the 
noise levels described in earlier in this section; the other was with the noise reduced to almost a 
negligible level in order to separate the effects of confusion. The image size for the distributed 
source simulation is 8192" x 8192". As mentioned in the previous section, we expect that the 
confusion is important, especially for the cases that 7 is greater than 2.5. In figure 6, we show 
an example of the simulated images with the normal level of noise. 

We carried out aperture photometry on the simulated images using SExtractor software 
v2.0.0 (Bertin, Arnouts 1996). Some influential parameters were optimised for better detection 
of the source, while the remaining were left intact as default values. We set the threshold in 
the source detection and the analysis as 3, the size of the photometric aperture as FWHM of 
beam pattern, and we did not apply a filter for detection. In order to calibrate the output flux, 
we used the five brightest input sources. 

Figure 7 shows the distribution of the Sont/Sin as a function of Sin, where Sin and Sont 
denote the input flux and the flux obtained by photometry. In the upper-left panel of figure 
7, we assumed that it is for the case with 7 = 1.5 and Nq = 10, and negligible contribution 
of photon and readout noise. We also assumed that detected source corresponds to the input 
source if the position of the detected source lies within 9" for SW bands and 15" for LW bands 
from the input source location. We found very good correlation between the input and output 
fluxes, and hence can conclude that the confusion noise is also negligible for this case. 

The noise added results for the case with 7 = 1.5 are shown in the upper- left panel of 
figure 8. The flux uncertainty becomes significant near the estimated detection limits due to 
photon and readout noise. Below the detection limit, most of the detected sources have an 
output flux greater than the input flux: This is simply because detection can be possible only 
when positive noises have been added to the source. 

The results with more crowded sources (i.e., 7 = 2.5 and 7 = 3.0) are shown in the lower 
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Figure06a. ps 



Figure06b . ps 



Fig. 6. Simulated images in tlic WIDE-L band for distributed sources. Wc generated the distributed 
sources according to the cases of 7 1.5 and Nf) ~ 10 (left panel) and 7 = 2.5 and A^q — 10 (right panel). 
Photon and readout noises are added in these images. 
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Figure07a. ps 



FigureOTb . ps 



FigureOTc . ps 



Figure07d.ps 



FigureOSa. ps 



FigureOSb . ps 



FigureOSc . ps 



Figure08d.ps 



panels of figure 7 for negligible noises, and figure 8 for normal noises. Even with negligible 
noises, we find that there are large deviations of the output fluxes from the input fluxes. Thus, 
the flux uncertainties are mostly caused by the source confusion shown in figure 7. Similar to 
the case dominated by the photon and readout noise shown in the upper-right panel of figure 
8, ^out is systematically overestimated for sources below the theoretical confusion limits. Such 
an upward bias was caused by source confusion; many of the detected sources contain fainter 
sources within the beam. Actually, the significant upward bias is partially due to the parameter, 
i.e., threshold, set in SExtrator. First, SExtractor estimates the background fluctuation from 
each local area. Because we reduce the noise below a negligible level, the calculated background 
fluctuations are mainly due to many dim sources. The detected sources at low flux surely have 
a flux above the fluctuation times the threshold; these detected sources cause a significant 
upward bias. In the case of heavy confusion, the trend of the boosted flux (see the lower panels 
of figure 8) is very similar to the case without noises (see the lower panels of figure 7), which 
means that the faint sources work as the dominant noise. 

Figure 9 shows the integrated source count results. For a comparison, we also plot the 
input source distribution. In the case of weak source confusion (i.e., 7 = 1.5 and A'^o = 10) 
(upper-left panel of figure 9), the source count from a simulated image follows the input source 
distribution well, except for the faint ends dominated by photon and readout noise. However, 
the lower panels of figure 9 show that the source distribution deviated from the input one due 
to source confusion. The location of the estimated confusion limit of table 4 is also shown in 
this figure. The observed slope is significantly different from the input slope. The output slope 
can be 1.5-times larger than the input slope in the case of a crowded source distribution. 

As we mentioned in subsubsection 4.2.2, we generated crowded fields for the case of 
7 = 1.5 by simply increasing Nq by a large factor, i.e., Ao(> 100 mJy) = 100 and the case of 
7 = 2.5 and A'o = 10, to exclude photon and readout noise in order to check the effect of pure 
source confusion. Because there are no significant difference between the case with and without 
the photon and readout noise, as can be seen in the upper-right panel of figures 7 and 8, we show 
the source count result with the photon and readout noise in the upper-right panel of figure 9 in 
comparison with the less-crowded case (upper left). Clearly, the confusion becomes important 
at around S = 100 mJy for WIDE-L according to a theoretical calculation, but the slope does 
not change. The change in the slope appears to occur only when the underlying A^(> S) varies 
rather steeply on S. The lower-right panel of figure 9 shows the case of 7 = 3.0 and Nq = 60, 
including the photon and readout noise. The slope of the source count is significantly changed 
by the heavy confusion, and the source detection mainly depends on source confusion. 

4-3.2. Detection limits from simulations 

It is not easy to define the detection limits from the simulated data. Since the detection 
becomes increasingly difficult for sources below the detection hmits, we first define the 'detection 
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Table 5. Detection limits for distributed point sources without photon and readout noise. 





7= 1.5 


7 


= 1.5 


7 = 2.5 


7 = 3.0 




N* = 10 


N* 


= 100 


N* = 10 


N* = 60 


Band 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


WIDE-L 


no confusion 




100 


58 


355 


N170 


no confusion 




105 


61 


390 


WIDE-S 


no confusion 




45 


31 


305 


N60 


no confusion 




40 


30 


278 



* Nq{> 100 mJy). Number per square degree. 

correctness' such that the ratio of the number of correctly detected sources to the number of 
detected sources from the photometry. We assume that the flux of the correctly detected source 
is the measured flux from the photometry, and agrees with the input flux within a 20% error. 
The detection correctness can be near unity for sources well beyond the detection limit, and 
goes down rapidly below the detection limit. We find that the detection correctness reaches 
around 0.7 at the estimated detection limit of a single scan. We thus define the location of the 
70% detection correctness as the detection limit in our simulated data. 

Figure 10 show a plot of the detection correctness with the photon and readout noise. 
We first attempted to estimate the detection limit purely due to source confusion. We arbi- 
trarily suppressed the photon and readout noise by a factor of 100 so that the noise-dominated 
detection limit would become much less than the lower limit of the source flux of 10 mJy. The 
resulting detection limits, estimated based on the detection correctness, are summarized in 
table 5. Under this condition, because the source detection is affected by the source confusion 
and the photometric accuracy, we could obtain similar detection limits in both narrow and wide 
bands. These numbers are similar to those in table 4, except for 7 = 1.5, where the detection 
correctness remains larger than 0.7, even for the faintest sources and for the case of the crowded 
source distribution. This means that the confusion is not important for 7 = 1.5 and A'^o(> 100 
mJy) = 10. 

Table 6 shows the estimates of combined detection limits where the readout noise, the 
photon noise, and the confusion noise are considered. Since the confusion is not important for 
the case of 7 = 1.5 and Nq = 10, the detection limit is purely determined by the photon and 
readout noise. For the case of 7 = 1.5 and A'o = 100, 7 = 2.5, and 7 = 3.0, both the source 
confusion and the other noises contribute to the detection limits. The combined detection limits 
for this case exceeds both the noise dominated result (table 2) and source confusion dominated 
result (table 4). In the case of 7 = 3.0 and Nq = 60, we cannot exactly determine the detection 
limits because the severe confusion makes the source detection difficult. Too many sources 
(i.e., 7 = 1.5 and Nq = 100, 7 = 3.0 and Nq = 60) also act as the large amount of the photon 
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Table 6. Detection limits for distributed point sources with photon and readout noise, taking account of the effects of 
the performance of the entire system, the brightness of the sky, the telescope emission, and the distribution of sources. 



Band 


7-1.5 

iVo* - 10 
(mJy) 


7-1.5 

iVo* = 100 
(mJy) 


7-2.5 

iVo* = 10 
(mJy) 


7-3.0 
- 60 
(mJy) 


WIDE-L 


26 


125 


68 


440 


N170 


66 


135 


115 


442 


WIDE-S 


21 


82 


40 


310 


N60 


49 


92 


63 


280 



* No{> lOOmJy). Number per square degree. 

noise, which affects in raising the detection hmit. Therefore, accurate photometry could be an 
additional important factor for approaching the theoretical confusion limit in these cases. 

Matsuhara et al. (2000) analysed the ISO data obtained for the high density case (7 = 3.0 
and A'^o = 60) from the fluctuation analysis method, which is different from our photometric 
method. Because they assumed that the fluctuation is mainly caused by unresolved faint point 
sources, they could count the number of sources, even in a low flux range. 

5. Summary 

We have written observing simulation software, 'FISVI', for an upcoming infrared survey 
mission, ASTRO-F. Utilizing this software, we have estimated the performance of the Far 
Infrared Surveyor (FIS) onboard ASTRO-F for ideal conditions. We can carry out scanning 
simulations with a reasonable amount of computing resources by introducing the Compiled 
PSF. The software can be used to generate virtual data sets for a data-reduction pipeline. 

We estimated the detection limits under various circumstances. For the case of a non- 
crowded source distribution, the readout noise is usually more important than the photon noise 
for dark patches of the sky by a factor of 1.3 to 2.5. This means that the bright parts of the 
sky can be easily dominated by photon noise. The emission from the telescope is less than 
the interstellar background as long as the telescope temperature remains less than 6 K, but it 
could contribute significantly to the long-wavelength band if the temperature becomes larger 
than 6.5 K (see figure 4). 

In crowded fields, source confusion becomes important in identifying sources. The detec- 
tion correctness becomes smaller for fainter sources. We have defined the confusion limit in such 
a way that the number of correctly detected sources within a 20% error becomes larger than 
70% of the number of detected sources from photometry. Such a definition of the confusion- 
dominated detection limit gives very similar values of the confusion limit based on a simple 
formula. The source confusion becomes larger than the detection limits by photon and read- 
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out noise only if the number of faint sources becomes much larger than a simple extension 
of the IRAS source counts down to around 10 mJy, assuming no luminosity or density evolu- 
tion. Recent models of source counts based on ISO and SCUBA observations (Matsuhara et 
al. 2000; Dole et al. 2001; Franceschini et al. 2001; Pearson 2001), however, predict the source 
distribution that is subject to significant confusion at the longest wavelength band (WIDE-L). 
Other bands appear to be noise-limited. The source confusion also could change the slope in 
logA^-logS" plots. 

In this paper, we have made many simplifying assumptions concerning the sky condi- 
tions. The actual sky brightness varies from place to place. The overall statistics of the galaxy 
counts should be significantly influenced by irregularities of the sky backgrounds. Also, in order 
to understand cosmo logical effects, we will consider various types of SED, the luminosity func- 
tion, and the redshift distribution. The current version of FISVI does not take into account 
more complicated behaviors of the detectors. These issues will be discussed in forthcoming 
papers. 
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Appendix 1. Compiled PSF 

A. 1.1. PSF Convolution 

The PSF of ASTRO-F/FIS, including the entire optical path, was computed using the 
ZEMAX optical simulation software package (Focus Software, Inc.). The resulting PSF at 
A = 200 /im is shown in figure 11, together with a circular aperture Airy pattern. The difference 
between the simulated PSF and the Airy pattern is very small, but noticeable. The simulated 
PSF is slightly narrower than the Airy pattern, and the side-lobe is more significant. Since 
FIS detectors do not lie on the optical axis of ASTRO-F, the PSF is slightly elongated with an 
ellipticity of ~ 0.05, but we assume the circular PSF in the present simulations. Since the FIS 
covers a wide range of wavelengths, the PSFs have been computed from 40 to 200 /im at 5 /xm 
intervals. 

Using the simulated PSF, we first obtain the PSF-convolved image Ix,i on the focal plane 
at wavelength A, contributed solely by the z-th point source: 

IxAr) = Fx,^hx{r■y^ (Al) 

and 
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Fig. 11. PSF of the ASTRO-F/FIS at 200 in a linear scale (left panel) and a logarithmic scale (right 
panel). The solid line shows the PSF simulated by using ZEMAX and the dotted line shows the Airy 
pattern with the assumption of a single circular aperture system. 
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hxir-y^)dQ, (A2) 

where r is the position vector on the focal plane, Fx^i is the flux density (at the wavelength A) 
of the 2-th source, and hx{r;r^j) is the simulated PSF at wavelength A located centered at the 
position of the i-th source r^. The PSF is normalised in such a way that the integration over 
the entire solid angle becomes unity. The intensity distribution on the focal plane, Ix{r), can 
then be obtained by 

Ix{r) = Y.Fx,hx{r;r'^). (A3) 

i 

A. 1.2. Filter Transmittance and Detector Response 

As the detector sweeps the sky, it integrates the charge generated by photons that fall 
onto the detector. For a given intensity distribution on the focal plane, Ix{f), the power, Px{f), 
at the wavelength interval dX is 

Px{r)dX= [ Ix{r) At,iT{X) dndX, (A4) 



where Atd is the effective collecting area of the telescope, and r(A) is the filter transmittance 
along the photon path within FIS (Takahashi et al. 2000). The integration is performed over 
the solid angle subtended by the pixel. 

The detector transforms the photons into charges. The total charge, D, integrated from 
ti to t2 is 

D{t,^t2) = f r Px{r{t))^{X) dtdX, (A5) 
JxJti 

where ^(A) is the detector response function in units of A W~^. We use the following convention: 

e(A)=eof(A), (A6) 

where is ^ constant in units of AW~^ and ^ is a function normalised to unity at the peak 
value for SW (short wavelength) and LW (long wavelength) detectors. 

The normalised detector response functions, C,, of LW and SW bands are shown in 
figures 12 and 13, respectively. We use these curves and the measured detector responsivity, C,^, 
to determine the normalisation constant, ,^o- Measurements are done using a blackbody source, 
a filter that cuts off the photons below a certain wavelength, a Winston cone, and a detector 
in a perfectly reflecting cavity. The LW detector has long wavelength cut-off at 200 /xm and 
SW detector at 110 fim. A low-pass filter was used to cut off the photons at wavelength below 
the FIS band. The short wavelength limits were 140 /im for the LW detector and 40 /im for 
the SW detector. The measured responsivity is represented by 

^ _^ IxaX)BxiT) dX 

IxBx{T)dX ' (^^^ 
where Bx{T) is the Planck function at the temperature T. In this estimation, we use T = 40 K. 
From the measured value of ~ 20 A W~^ for LW, and ~ 7 A W~^ for SW, we can determine 
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Fig. 12. Filter transmission, t{X), (dashed lines) and the detector's response function, ^(A), (dotted lines) 
for the N60 band (left) and the Wide-S band (right). The combined responsivities are shown as solid lines 
in arbitrary units. 

the normalisation constant, ^o- The normahsation constants are = 30 A for the LW and 
^0 = 10 A for the SW detectors, respectively. 

A. 1.3. Compiled PSF 

If we use the same spectral energy distribution (SED) for each source, the flux density 
of the source can be defined as 



(A8) 



where 5'a is the spectral energy distribution (SED) normalised to unity over the wavelength 
band and Ti is the flux integrated over the bandwidth. We can rewrite equation (A3) as 
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Fig. 13. Same as figure 12, exeept for tire N170 band (left) and tire WIDE-L band (right). 
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(A9) 



Since A is independent of r and r^, we can introduce a new function, H{r\ r'J, by integrating 
over the wavelength as 



We define this H[r\r^^ as the 'Compiled PSF'. If we perform convolution to the image plane 
by using this Compiled PSF, we can avoid repeated wavelength integration. Finally, equation 
(A5) can be rewritten as 



This concept of the Compiled PSF is effective only when the number of SED type is limited. 
The calculation time is reduced by a factor of Nx by using the Compiled PSF, where Nx is 
the number of wavelength grids. With a wavelength interval of AA = 5 /im, a typical A^a lies 
between 10 and 20. In order to carry out simulations over four square degrees in the WIDE- 
S band, we need about 15 hours of computing time with Pentium IV 1 GHz machines. By 
introducing Compiled PSF, we can accomplish such a simulation within an hour. 

A. 1.4- Spectral Energy Distribution of the Sources 

We expect that the majority of faint point sources detected by the ASTRO-F/FIS will 
be external galaxies. Each object will have its own SED, but most extragalactic point sources in 
the infrared band can be classified into four types of galaxies, i.e., the cirrus type representing 
typical spiral galaxies, the M 82 type starbursts, the Arp 220 type starbursts and the AGN 
dust torus type (Rowan- Robinson 2001). Four Compiled PSFs are required to accommodate 
these four types of SEDs in the simulations. The observed SEDs are further affected by the 
redshifts. We need redshifted-dependent SEDs for each type of source. 

We expect that the Compiled PSF will be changed with the SED types and the redshift 
for wide bands, but the difference was found to be very small, even for the WIDE-S and WIDE- 
L bands, as shown in figure 14. Since our main purpose is to examine the general performance 
of the ASTRO-F/FIS, we concentrate on simple models for the nature of the sources. We will 
deal with the SED types of sources, redshift distributions, and the luminosity function in the 
next paper in order to understand the cosmological model and the galaxy evolution through 
the observing simulation. Though the difference between the Compiled PSFs computed from 
the flat SED and other SEDs is severe at some extreme cases (~ 10% difference over the area), 
we use the Compiled PSF computed for galaxies with the fiat SED in the present paper (i.e.. 
Fx = constant) (see figure 14). In the fiat SED's case, the Compiled PSF does not depend on 
the redshift. 




(AlO) 




(All) 
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Fig. 14. Normalised Compiled PSFs in the WIDE-L band. The left panel is the Compiled PSFs over 
SED for redshift 1.0 and the right panel is the Compiled PSFs over redshift for the cirrus type. For a 
comparison, we also plot the Compiled PSF computed for galaxies with the flat SED used in this work. 
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Appendix 2. Procedures of Scanning and Data Sampling 



A PSF-convolved image is generated on grids wliere tlie scanning procedure is performed. 
To scan a PSF-convolved image, we need to know the position of the detector pixels. We set the 
array of the starting point to scan on the x (cross-scan direction) and the y (in-scan direction) 
frame in the image. The FIS detector arrays have 2 or 3 rows and 15 or 20 columns, and is 
tilted by an angle 6 = 26. °5 from the cross-scan direction in order to assure Nyquist sampling 
(Takahashi et al. 2000; Matsuura et al. 2001). We denote i as the index for the sampling 
sequence, and j and k as the indices for the row and column of the detector array, respectively 
(see figure 15). By denoting (xq, yo) as the position vector of the center of upper left pixel of 
the array at the beginning of the scan (i.e., i = j = k = 0), we have the following formulae for 
the position vectors of the {j,k) pixel at the {i + l)-th sampling: 

x{i,j, k) = xo + p{k cos9 + j sin 9) (A12) 

and 

y{i,j, /c) = yo + ivAt + p{j cos 9 — k sin 9) , (A13) 

where p is the size of the pixel pitch (see table 1), f is the scanning angular speed (which is 
3.60 arcmin s"^) of the satellite, and At is the increment of the detector motion in the scan 
direction during the sampling interval. Note that the x position of each pixel does not depend 
on i in this coordinate system. We show one example for the passage of the detector in figure 
16. 

Appendix 3. Integrating over the Detector Pixel 

The integration of equation (All) over r2pixci was carried out by summing up the image 
convolved with Compiled PSF on fine grids. The image convolved with Compiled PSF was 
constructed on grids of 4" resolution, but the accuracy of the fipixei integration was not good 
enough on such grids (~ a few percent error), partly because of the tilted configuration of the 
detector arrays. In order to improve the accuracy of the integration, we laid finer grids over the 
area where the integration would be performed. We were able to reduce the integration error 
down to 1% by taking a three-times finer grid over the integration area. If we use a smaller 
grid, we can improve the accuracy of the fiux and the position further, but we would need more 
computing time. 

The time integration of equation (All) was made by dividing one sampling interval to 
shorter subsampling intervals in order to mimic the continuous scanning of the detector and 
applying the trapezoidal rule to the subsampled time series data. As the detector moves, one 
detector pixel integrates the signal during the subsample interval (see figure 17). The number 
of subsample determines the resolution of integrated signal values. The sampling rate of 15.2 
Hz for LW bands corresponds to 14f'2 which is much smaller than the pixel size, and we found 
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Fig. 15. Layout of the detector array for the N170 band and definition of the scan directions. 
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Fig. 16. Passage of the detector for the WIDE-L band. We display the footprints of the detector pixels 
scanned three times. 
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Figure 17 .ps 



Fig. 17. Each sampling is composed of several subsamples in order to ensure accurate integration over 
the region where the intensity varies. sl-s5 mean the subsamples. In actual simulations, we used only two 
subsamples. 

that we need only two subsamples to ensure the integration accuracy over time becomes smaller 
than 1%. 
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